Quantitative Assessment of Tumor Contact with Neurogenic Zones and Its Effects on Survival: Insights beyond Traditional Predictors.

So far, the cellular origin of glioblastoma (GBM) needs to be determined, with prevalent theories suggesting emergence from transformed endogenous stem cells. Adult neurogenesis primarily occurs in two brain regions: the subventricular zone (SVZ) and the subgranular zone (SGZ) of the hippocampal dentate gyrus. Whether the proximity of GBM to these neurogenic niches affects patient outcome remains uncertain. Previous studies often rely on subjective assessments, limiting the reliability of those results. In this study, we assessed the impact of GBM's relationship with the cortex, SVZ and SGZ on clinical variables using fully automated segmentation methods. In 177 glioblastoma patients, we calculated optimal cutpoints of minimal distances to the SVZ and SGZ to distinguish poor from favorable survival. The impact of tumor contact with neurogenic zones on clinical parameters, such as overall survival, multifocality, MGMT promotor methylation, Ki-67 and KPS score was also examined by multivariable regression analysis, chi-square test and Mann-Whitney-U. The analysis confirmed shorter survival in tumors contacting the SVZ with an optimal cutpoint of 14 mm distance to the SVZ, separating poor from more favorable survival. In contrast, tumor contact with the SGZ did not negatively affect survival. We did not find significant correlations with multifocality or MGMT promotor methylation in tumors contacting the SVZ, as previous studies discussed. These findings suggest that the spatial relationship between GBM and neurogenic niches needs to be assessed differently. Objective measurements disprove prior assumptions, warranting further research on this topic.

LST-AI: A deep learning ensemble for accurate MS lesion segmentation.

Automated segmentation of brain white matter lesions is crucial for both clinical assessment and scientific research in multiple sclerosis (MS). Over a decade ago, we introduced an engineered lesion segmentation tool, LST. While recent lesion segmentation approaches have leveraged artificial intelligence (AI), they often remain proprietary and difficult to adopt. As an open-source tool, we present LST-AI, an advanced deep learning-based extension of LST that consists of an ensemble of three 3D U-Nets. LST-AI explicitly addresses the imbalance between white matter (WM) lesions and non-lesioned WM. It employs a composite loss function incorporating binary cross-entropy and Tversky loss to improve segmentation of the highly heterogeneous MS lesions. We train the network ensemble on 491 MS pairs of T1-weighted and FLAIR images, collected in-house from a 3T MRI scanner, and expert neuroradiologists manually segmented the utilized lesion maps for training. LST-AI also includes a lesion location annotation tool, labeling lesions as periventricular, infratentorial, and juxtacortical according to the 2017 McDonald criteria, and, additionally, as subcortical. We conduct evaluations on 103 test cases consisting of publicly available data using the Anima segmentation validation tools and compare LST-AI with several publicly available lesion segmentation models. Our empirical analysis shows that LST-AI achieves superior performance compared to existing methods. Its Dice and F1 scores exceeded 0.62, outperforming LST, SAMSEG (Sequence Adaptive Multimodal SEGmentation), and the popular nnUNet framework, which all scored below 0.56. Notably, LST-AI demonstrated exceptional performance on the MSSEG-1 challenge dataset, an international WM lesion segmentation challenge, with a Dice score of 0.65 and an F1 score of 0.63-surpassing all other competing models at the time of the challenge. With increasing lesion volume, the lesion detection rate rapidly increased with a detection rate of >75% for lesions with a volume between 10 mm3 and 100 mm3. Given its higher segmentation performance, we recommend that research groups currently using LST transition to LST-AI. To facilitate broad adoption, we are releasing LST-AI as an open-source model, available as a command-line tool, dockerized container, or Python script, enabling diverse applications across multiple platforms.

Single brain metastases - prognostic factors and impact of residual tumor burden on overall survival.

Background: Brain metastases (BM) are a common and challenging issue, with their incidence on the rise due to advancements in systemic therapies and increased patient survival. Most patients present with single BM, some of them without any further extracranial metastasis (i.e., solitary BM). The significance of postoperative intracranial tumor volume in the treatment of singular and solitary BM is still debated. Objective: This study aimed to determine the impact of resection and postoperative tumor burden on overall survival (OS) in patients with single BM. Methods: Patients with surgically treated single BM between 04/2007-01/2020 were retrospectively included. Residual tumor burden (RTB) was determined by manual segmentation of early postoperative brain MRI (72 h). Survival analyses were performed using Kaplan-Meier estimates for univariate analysis and Cox regression proportional hazards model for multivariate analysis, using preoperative Karnofsky performance status scale (KPSS), age, sex, RTB, incomplete resection and singular/solitary BM as covariates. Results: 340 patients were included, median age 64 years (54-71). 119 patients (35%) had solitary BM, 221 (65%) singular BM. Complete resection (RTB=0) was achieved in 73%, median preoperative tumor burden was 11.2 cm3 (5-25), and RTB 0 cm3 (0-0.2). Median OS of patients with singular BM was 13 months (4-33) vs 20 months (5-92) for solitary BM; p=0.062. Multivariate analysis revealed singular BM as independent risk factor for poorer OS: HR 1.840 (1.202-2.817), p=0.005. Complete vs. incomplete resection showed no significant OS difference (13 vs. 13 months, p=0.737). When focusing on solitary BM, complete resection led to a longer OS than incomplete resection (21 vs. 8 months), without statistical significance(p=0.250). Achieving RTB=0 resulted in higher OS for patients with solitary BM compared to singular BM (21 vs. 12 months, p=0.027). Patients who received postoperative radiotherapy (RT) had significantly longer OS compared to those without it (14 vs. 4 months, p<0.001), with favorable OS in those receiving stereotactic radiosurgery (SRS) (15 months (3-42), p<0.001) or hypofractionated stereotactic radiotherapy (HSRT). Conclusion: When complete intracranial tumor resection RTB=0 is achieved, patients with solitary BM have a favorable outcome compared to singular BM. Singular BM was confirmed as independent risk factor. There is a strong presumption that complete resection leads to an improved oncological prognosis. Patients with solitary BM tend to benefit with a favorable outcome following complete resection. Hence, surgical resection should be considered as a treatment option for patients presenting with either no or minimal extracranial disease. Furthermore, the highly favorable impact of postoperative RT on OS was demonstrated and confirmed, especially with SRS or HSRT.

Quantitative susceptibility mapping in multiple sclerosis: A systematic review and meta-analysis.

BACKGROUND: Quantitative susceptibility mapping (QSM) is a quantitative measure based on magnetic resonance imaging sensitive to iron and myelin content. This makes QSM a promising non-invasive tool for multiple sclerosis (MS) in research and clinical practice. OBJECTIVE: We performed a systematic review and meta-analysis on the use of QSM in MS. METHODS: Our review was prospectively registered on PROSPERO (CRD42022309563). We searched five databases for studies published between inception and 30th April 2023. We identified 83 English peer-reviewed studies that applied QSM images on MS cohorts. Fifty-five included studies had at least one of the following outcome measures: deep grey matter QSM values in MS, either compared to healthy controls (HC) (k = 13) or correlated with the score on the Expanded Disability Status Scale (EDSS) (k = 7), QSM lesion characteristics (k = 22) and their clinical correlates (k = 17), longitudinal correlates (k = 11), histological correlates (k = 7), or correlates with other imaging techniques (k = 12). Two meta-analyses on deep grey matter (DGM) susceptibility data were performed, while the remaining findings could only be analyzed descriptively. RESULTS: After outlier removal, meta-analyses demonstrated a significant increase in the basal ganglia susceptibility (QSM values) in MS compared to HC, caudate (k = 9, standardized mean difference (SDM) = 0.54, 95 % CI = 0.39-0.70, I2 = 46 %), putamen (k = 9, SDM = 0.38, 95 % CI = 0.19-0.57, I2 = 59 %), and globus pallidus (k = 9, SDM = 0.48, 95 % CI = 0.28-0.67, I2 = 60 %), whereas thalamic QSM values exhibited a significant reduction (k = 12, SDM = -0.39, 95 % CI = -0.66--0.12, I2 = 84 %); these susceptibility differences in MS were independent of age. Further, putamen QSM values positively correlated with EDSS (k = 4, r = 0.36, 95 % CI = 0.16-0.53, I2 = 0 %). Regarding rim lesions, four out of seven studies, representing 73 % of all patients, reported rim lesions to be associated with more severe disability. Moreover, lesion development from initial detection to the inactive stage is paralleled by increasing, plateauing (after about two years), and gradually decreasing QSM values, respectively. Only one longitudinal study provided clinical outcome measures and found no association. Histological data suggest iron content to be the primary source of QSM values in DGM and at the edges of rim lesions; further, when also considering data from myelin water imaging, the decrease of myelin is likely to drive the increase of QSM values within WM lesions. CONCLUSIONS: We could provide meta-analytic evidence for DGM susceptibility changes in MS compared to HC; basal ganglia susceptibility is increased and, in the putamen, associated with disability, while thalamic susceptibility is decreased. Beyond these findings, further investigations are necessary to establish the role of QSM in MS for research or even clinical routine.

Noninvasive- and invasive mapping reveals similar language network centralities - A function-based connectome analysis.

BACKGROUND: Former comparisons between direct cortical stimulation (DCS) and navigated transcranial magnetic stimulation (nTMS) only focused on cortical mapping. While both can be combined with diffusion tensor imaging, their differences in the visualization of subcortical and even network levels remain unclear. Network centrality is an essential parameter in network analysis to measure the importance of nodes identified by mapping. Those include Degree centrality, Eigenvector centrality, Closeness centrality, Betweenness centrality, and PageRank centrality. While DCS and nTMS have repeatedly been compared on the cortical level, the underlying network identified by both has not been investigated yet. METHOD: 27 patients with brain lesions necessitating preoperative nTMS and intraoperative DCS language mapping during awake craniotomy were enrolled. Function-based connectome analysis was performed based on the cortical nodes obtained through the two mapping methods, and language-related network centralities were compared. RESULTS: Compared with DCS language mapping, the positive predictive value of cortical nTMS language mapping is 74.1%, with good consistency of tractography for the arcuate fascicle and superior longitudinal fascicle. Moreover, network centralities did not differ between the two mapping methods. However, ventral stream tracts can be better traced based on nTMS mappings, demonstrating its strengths in acquiring language-related networks. In addition, it showed lower centralities than other brain areas, with decentralization as an indicator of language function loss. CONCLUSION: This study deepens the understanding of language-related functional anatomy and proves that non-invasive mapping-based network analysis is comparable to the language network identified via invasive cortical mapping.

Human-Level Differentiation of Medulloblastoma from Pilocytic Astrocytoma: A Real-World Multicenter Pilot Study.

Medulloblastoma and pilocytic astrocytoma are the two most common pediatric brain tumors with overlapping imaging features. In this proof-of-concept study, we investigated using a deep learning classifier trained on a multicenter data set to differentiate these tumor types. We developed a patch-based 3D-DenseNet classifier, utilizing automated tumor segmentation. Given the heterogeneity of imaging data (and available sequences), we used all individually available preoperative imaging sequences to make the model robust to varying input. We compared the classifier to diagnostic assessments by five readers with varying experience in pediatric brain tumors. Overall, we included 195 preoperative MRIs from children with medulloblastoma (n = 69) or pilocytic astrocytoma (n = 126) across six university hospitals. In the 64-patient test set, the DenseNet classifier achieved a high AUC of 0.986, correctly predicting 62/64 (97%) diagnoses. It misclassified one case of each tumor type. Human reader accuracy ranged from 100% (expert neuroradiologist) to 80% (resident). The classifier performed significantly better than relatively inexperienced readers (p < 0.05) and was on par with pediatric neuro-oncology experts. Our proof-of-concept study demonstrates a deep learning model based on automated tumor segmentation that can reliably preoperatively differentiate between medulloblastoma and pilocytic astrocytoma, even in heterogeneous data.

Virtual reality-empowered deep-learning analysis of brain cells.

Automated detection of specific cells in three-dimensional datasets such as whole-brain light-sheet image stacks is challenging. Here, we present DELiVR, a virtual reality-trained deep-learning pipeline for detecting c-Fos+ cells as markers for neuronal activity in cleared mouse brains. Virtual reality annotation substantially accelerated training data generation, enabling DELiVR to outperform state-of-the-art cell-segmenting approaches. Our pipeline is available in a user-friendly Docker container that runs with a standalone Fiji plugin. DELiVR features a comprehensive toolkit for data visualization and can be customized to other cell types of interest, as we did here for microglia somata, using Fiji for dataset-specific training. We applied DELiVR to investigate cancer-related brain activity, unveiling an activation pattern that distinguishes weight-stable cancer from cancers associated with weight loss. Overall, DELiVR is a robust deep-learning tool that does not require advanced coding skills to analyze whole-brain imaging data in health and disease.

Intensity scaling of conventional brain magnetic resonance images avoiding cerebral reference regions: A systematic review.

BACKGROUND: Conventional brain magnetic resonance imaging (MRI) produces image intensities that have an arbitrary scale, hampering quantification. Intensity scaling aims to overcome this shortfall. As neurodegenerative and inflammatory disorders may affect all brain compartments, reference regions within the brain may be misleading. Here we summarize approaches for intensity scaling of conventional T1-weighted (w) and T2w brain MRI avoiding reference regions within the brain. METHODS: Literature was searched in the databases of Scopus, PubMed, and Web of Science. We included only studies that avoided reference regions within the brain for intensity scaling and provided validating evidence, which we divided into four categories: 1) comparative variance reduction, 2) comparative correlation with clinical parameters, 3) relation to quantitative imaging, or 4) relation to histology. RESULTS: Of the 3825 studies screened, 24 fulfilled the inclusion criteria. Three studies used scaled T1w images, 2 scaled T2w images, and 21 T1w/T2w-ratio calculation (with double counts). A robust reduction in variance was reported. Twenty studies investigated the relation of scaled intensities to different types of quantitative imaging. Statistically significant correlations with clinical or demographic data were reported in 8 studies. Four studies reporting the relation to histology gave no clear picture of the main signal driver of conventional T1w and T2w MRI sequences. CONCLUSIONS: T1w/T2w-ratio calculation was applied most often. Variance reduction and correlations with other measures suggest a biologically meaningful signal harmonization. However, there are open methodological questions and uncertainty on its biological underpinning. Validation evidence on other scaling methods is even sparser.

Impact of function-guided glioma treatment on oncological outcome in the elderly.

Introduction: Many patients with high-grade gliomas (HGG) are of older age. Research question: We hypothesize that pre- and intraoperative mapping and monitoring preserve functional status in elderly patients while gross total resection (GTR) is the aim, resulting in overall survival (OS) rates comparable to the general population with HGG. Material and methods: We subdivided a prospective cohort of 168 patients above 65 years with eloquent high-grade gliomas into four groups ([years/cases] 1: 65-69/58; 2: 70-74/47; 3: 75-79/43; 4: >79/20). All patients underwent preoperative noninvasive mapping, which was also used for decision-making, intraoperative neuromonitoring in 138 cases, direct cortical and/or subcortical motor mapping in 66 and 50 cases, and awake language mapping in 11 cases. Results: GTR and subtotal resection (STR) could be achieved in 65% and 28%, respectively. Stereotactic biopsy was performed in 8% of cases. Postoperatively, we found transient and permanent functional deficits in 13% and 11% of cases. Postoperative Karnofsky Performance Scale (KPS) did not differ between subgroups. Patients with long-term follow-up (51%) had a progression-free survival of 5.5 (1-47) months and an overall survival of 10.5 (0-86) months. Discussion and conclusion: The interdisciplinary glioma treatment in the elderly is less age-dependent but must be adjusted to the functional status. Function-guided surgical resections could be performed as usual, with maximal tumor resection being the primary goal. However, less network capacity in the elderly to compensate for deficits might cause higher rates of permanent deficits in this group of patients with more fast-growing malignant gliomas.

A Learnable Prior Improves Inverse Tumor Growth Modeling.

Biophysical modeling, particularly involving partial differential equations (PDEs), offers significant potential for tailoring disease treatment protocols to individual patients. However, the inverse problem-solving aspect of these models presents a substantial challenge, either due to the high computational requirements of model-based approaches or the limited robustness of deep learning (DL) methods. We propose a novel framework that leverages the unique strengths of both approaches in a synergistic manner. Our method incorporates a DL ensemble for initial parameter estimation, facilitating efficient downstream evolutionary sampling initialized with this DL-based prior. We showcase the effectiveness of integrating a rapid deep-learning algorithm with a high-precision evolution strategy in estimating brain tumor cell concentrations from magnetic resonance images. The DL-Prior plays a pivotal role, significantly constraining the effective sampling-parameter space. This reduction results in a fivefold convergence acceleration and a Dice-score of 95.

Deep Learning to Differentiate Benign and Malignant Vertebral Fractures at Multidetector CT.

Background Differentiating between benign and malignant vertebral fractures poses diagnostic challenges. Purpose To investigate the reliability of CT-based deep learning models to differentiate between benign and malignant vertebral fractures. Materials and Methods CT scans acquired in patients with benign or malignant vertebral fractures from June 2005 to December 2022 at two university hospitals were retrospectively identified based on a composite reference standard that included histopathologic and radiologic information. An internal test set was randomly selected, and an external test set was obtained from an additional hospital. Models used a three-dimensional U-Net encoder-classifier architecture and applied data augmentation during training. Performance was evaluated using the area under the receiver operating characteristic curve (AUC) and compared with that of two residents and one fellowship-trained radiologist using the DeLong test. Results The training set included 381 patients (mean age, 69.9 years ± 11.4 [SD]; 193 male) with 1307 vertebrae (378 benign fractures, 447 malignant fractures, 482 malignant lesions). Internal and external test sets included 86 (mean age, 66.9 years ± 12; 45 male) and 65 (mean age, 68.8 years ± 12.5; 39 female) patients, respectively. The better-performing model of two training approaches achieved AUCs of 0.85 (95% CI: 0.77, 0.92) in the internal and 0.75 (95% CI: 0.64, 0.85) in the external test sets. Including an uncertainty category further improved performance to AUCs of 0.91 (95% CI: 0.83, 0.97) in the internal test set and 0.76 (95% CI: 0.64, 0.88) in the external test set. The AUC values of residents were lower than that of the best-performing model in the internal test set (AUC, 0.69 [95% CI: 0.59, 0.78] and 0.71 [95% CI: 0.61, 0.80]) and external test set (AUC, 0.70 [95% CI: 0.58, 0.80] and 0.71 [95% CI: 0.60, 0.82]), with significant differences only for the internal test set (P < .001). The AUCs of the fellowship-trained radiologist were similar to those of the best-performing model (internal test set, 0.86 [95% CI: 0.78, 0.93; P = .39]; external test set, 0.71 [95% CI: 0.60, 0.82; P = .46]). Conclusion Developed models showed a high discriminatory power to differentiate between benign and malignant vertebral fractures, surpassing or matching the performance of radiology residents and matching that of a fellowship-trained radiologist. © RSNA, 2024 See also the editorial by Booz and D'Angelo in this issue.

Toward image-based personalization of glioblastoma therapy: A clinical and biological validation study of a novel, deep learning-driven tumor growth model.

Background: The diffuse growth pattern of glioblastoma is one of the main challenges for accurate treatment. Computational tumor growth modeling has emerged as a promising tool to guide personalized therapy. Here, we performed clinical and biological validation of a novel growth model, aiming to close the gap between the experimental state and clinical implementation. Methods: One hundred and twenty-four patients from The Cancer Genome Archive (TCGA) and 397 patients from the UCSF Glioma Dataset were assessed for significant correlations between clinical data, genetic pathway activation maps (generated with PARADIGM; TCGA only), and infiltration (Dw) as well as proliferation (ρ) parameters stemming from a Fisher-Kolmogorov growth model. To further evaluate clinical potential, we performed the same growth modeling on preoperative magnetic resonance imaging data from 30 patients of our institution and compared model-derived tumor volume and recurrence coverage with standard radiotherapy plans. Results: The parameter ratio Dw/ρ (P < .05 in TCGA) as well as the simulated tumor volume (P < .05 in TCGA/UCSF) were significantly inversely correlated with overall survival. Interestingly, we found a significant correlation between 11 proliferation pathways and the estimated proliferation parameter. Depending on the cutoff value for tumor cell density, we observed a significant improvement in recurrence coverage without significantly increased radiation volume utilizing model-derived target volumes instead of standard radiation plans. Conclusions: Identifying a significant correlation between computed growth parameters and clinical and biological data, we highlight the potential of tumor growth modeling for individualized therapy of glioblastoma. This might improve the accuracy of radiation planning in the near future.

(Predictable) performance bias in unsupervised anomaly detection.

BACKGROUND: With the ever-increasing amount of medical imaging data, the demand for algorithms to assist clinicians has amplified. Unsupervised anomaly detection (UAD) models promise to aid in the crucial first step of disease detection. While previous studies have thoroughly explored fairness in supervised models in healthcare, for UAD, this has so far been unexplored. METHODS: In this study, we evaluated how dataset composition regarding subgroups manifests in disparate performance of UAD models along multiple protected variables on three large-scale publicly available chest X-ray datasets. Our experiments were validated using two state-of-the-art UAD models for medical images. Finally, we introduced subgroup-AUROC (sAUROC), which aids in quantifying fairness in machine learning. FINDINGS: Our experiments revealed empirical "fairness laws" (similar to "scaling laws" for Transformers) for training-dataset composition: Linear relationships between anomaly detection performance within a subpopulation and its representation in the training data. Our study further revealed performance disparities, even in the case of balanced training data, and compound effects that exacerbate the drop in performance for subjects associated with multiple adversely affected groups. INTERPRETATION: Our study quantified the disparate performance of UAD models against certain demographic subgroups. Importantly, we showed that this unfairness cannot be mitigated by balanced representation alone. Instead, the representation of some subgroups seems harder to learn by UAD models than that of others. The empirical "fairness laws" discovered in our study make disparate performance in UAD models easier to estimate and aid in determining the most desirable dataset composition. FUNDING: European Research Council Deep4MI.

Amino acid PET vs. RANO MRI for prediction of overall survival in patients with recurrent high grade glioma under bevacizumab therapy.

PURPOSE: To summarize evidence on the comparative value of amino acid (AA) PET and conventional MRI for prediction of overall survival (OS) in patients with recurrent high grade glioma (rHGG) under bevacizumab therapy. METHODS: Medical databases were screened for studies with individual data on OS, follow-up MRI, and PET findings in the same patient. MRI images were assessed according to the RANO criteria. A receiver operating characteristic curve analysis was used to predict OS at 9 months. RESULTS: Five studies with a total of 72 patients were included. Median OS was significantly lower in the PET-positive than in the PET-negative group. PET findings predicted OS with a pooled sensitivity and specificity of 76% and 71%, respectively. Corresponding values for MRI were 32% and 82%. Area under the curve and sensitivity were significantly higher for PET than for MRI. CONCLUSION: For monitoring of patients with rHGG under bevacizumab therapy, AA-PET should be preferred over RANO MRI.

Navigating Post-Operative Outcomes: A Comprehensive Reframing of an Original Graded Prognostic Assessment in Patients with Brain Metastases.

BACKGROUND: Graded Prognostic Assessment (GPA) has been proposed for various brain metastases (BMs) tailored to the primary histology and molecular profiles. However, it does not consider whether patients have been operated on or not and does not include surgical outcomes as prognostic factors. The residual tumor burden (RTB) is a strong predictor of overall survival. We validated the GPA score and introduced "volumetric GPA" in the largest cohort of operated patients and further explored the role of RTB as an additional prognostic factor. METHODS: A total of 630 patients with BMs between 2007 and 2020 were included. The four GPA components were analyzed. The validity of the original score was assessed using Cox regression, and a modified index incorporating RTB was developed by comparing the accuracy, sensitivity, specificity, F1-score, and AUC parameters. RESULTS: GPA categories showed an association with survival: age (p < 0.001, hazard ratio (HR) 2.9, 95% confidence interval (CI) 2.5-3.3), Karnofsky performance status (KPS) (p < 0.001, HR 1.3, 95% CI 1.2-1.5), number of BMs (p = 0.019, HR 1.4, 95% CI 1.1-1.8), and the presence of extracranial manifestation (p < 0.001, HR 3, 95% CI 1.6-2.5). The median survival for GPA 0-1 was 4 months; for GPA 1.5-2, it was 12 months; for GPA 2.5-3, it was 21 months; and for GPA 3.5-4, it was 38 months (p < 0.001). RTB was identified as an independent prognostic factor. A cut-off of 2 cm3 was used for further analysis, which showed a median survival of 6 months (95% CI 4-8) vs. 13 months (95% CI 11-14, p < 0.001) for patients with RTB > 2 cm3 and <2 cm3, respectively. RTB was added as an additional component for a modified volumetric GPA score. The survival rates with the modified GPA score were: GPA 0-1: 4 months, GPA 1.5-2: 7 months, GPA 2.5-3: 18 months, and GPA 3.5-4: 34 months. Both scores showed good stratification, with the new score showed a trend towards better discrimination in patients with more favorable prognoses. CONCLUSION: The prognostic value of the original GPA was confirmed in our cohort of patients who underwent surgery for BM. The RTB was identified as a parameter of high prognostic significance and was incorporated into an updated "volumetric GPA". This score provides a novel tool for prognosis and clinical decision making in patients undergoing surgery. This method may be useful for stratification and patient selection for further treatment and in future clinical trials.

Personalized Predictions of Glioblastoma Infiltration: Mathematical Models, Physics-Informed Neural Networks and Multimodal Scans.

Predicting the infiltration of Glioblastoma (GBM) from medical MRI scans is crucial for understanding tumor growth dynamics and designing personalized radiotherapy treatment plans.Mathematical models of GBM growth can complement the data in the prediction of spatial distributions of tumor cells. However, this requires estimating patient-specific parameters of the model from clinical data, which is a challenging inverse problem due to limited temporal data and the limited time between imaging and diagnosis. This work proposes a method that uses Physics-Informed Neural Networks (PINNs) to estimate patient-specific parameters of a reaction-diffusion PDE model of GBM growth from a single 3D structural MRI snapshot. PINNs embed both the data and the PDE into a loss function, thus integrating theory and data. Key innovations include the identification and estimation of characteristic non-dimensional parameters, a pre-training step that utilizes the non-dimensional parameters and a fine-tuning step to determine the patient specific parameters. Additionally, the diffuse domain method is employed to handle the complex brain geometry within the PINN framework. Our method is validated both on synthetic and patient datasets, and shows promise for real-time parametric inference in the clinical setting for personalized GBM treatment.

LST-AI: a Deep Learning Ensemble for Accurate MS Lesion Segmentation.

Automated segmentation of brain white matter lesions is crucial for both clinical assessment and scientific research in multiple sclerosis (MS). Over a decade ago, we introduced an engineered lesion segmentation tool, LST. While recent lesion segmentation approaches have leveraged artificial intelligence (AI), they often remain proprietary and difficult to adopt. As an open-source tool, we present LST-AI, an advanced deep learning-based extension of LST that consists of an ensemble of three 3D-UNets. LST-AI explicitly addresses the imbalance between white matter (WM) lesions and non-lesioned WM. It employs a composite loss function incorporating binary cross-entropy and Tversky loss to improve segmentation of the highly heterogeneous MS lesions. We train the network ensemble on 491 MS pairs of T1w and FLAIR images, collected in-house from a 3T MRI scanner, and expert neuroradiologists manually segmented the utilized lesion maps for training. LST-AI additionally includes a lesion location annotation tool, labeling lesion location according to the 2017 McDonald criteria (periventricular, infratentorial, juxtacortical, subcortical). We conduct evaluations on 103 test cases consisting of publicly available data using the Anima segmentation validation tools and compare LST-AI with several publicly available lesion segmentation models. Our empirical analysis shows that LST-AI achieves superior performance compared to existing methods. Its Dice and F1 scores exceeded 0.62, outperforming LST, SAMSEG (Sequence Adaptive Multimodal SEGmentation), and the popular nnUNet framework, which all scored below 0.56. Notably, LST-AI demonstrated exceptional performance on the MSSEG-1 challenge dataset, an international WM lesion segmentation challenge, with a Dice score of 0.65 and an F1 score of 0.63-surpassing all other competing models at the time of the challenge. With increasing lesion volume, the lesion detection rate rapidly increased with a detection rate of >75% for lesions with a volume between 10mm3 and 100mm3. Given its higher segmentation performance, we recommend that research groups currently using LST transition to LST-AI. To facilitate broad adoption, we are releasing LST-AI as an open-source model, available as a command-line tool, dockerized container, or Python script, enabling diverse applications across multiple platforms.

Where is VALDO? VAscular Lesions Detection and segmentatiOn challenge at MICCAI 2021.

Imaging markers of cerebral small vessel disease provide valuable information on brain health, but their manual assessment is time-consuming and hampered by substantial intra- and interrater variability. Automated rating may benefit biomedical research, as well as clinical assessment, but diagnostic reliability of existing algorithms is unknown. Here, we present the results of the VAscular Lesions DetectiOn and Segmentation (Where is VALDO?) challenge that was run as a satellite event at the international conference on Medical Image Computing and Computer Aided Intervention (MICCAI) 2021. This challenge aimed to promote the development of methods for automated detection and segmentation of small and sparse imaging markers of cerebral small vessel disease, namely enlarged perivascular spaces (EPVS) (Task 1), cerebral microbleeds (Task 2) and lacunes of presumed vascular origin (Task 3) while leveraging weak and noisy labels. Overall, 12 teams participated in the challenge proposing solutions for one or more tasks (4 for Task 1-EPVS, 9 for Task 2-Microbleeds and 6 for Task 3-Lacunes). Multi-cohort data was used in both training and evaluation. Results showed a large variability in performance both across teams and across tasks, with promising results notably for Task 1-EPVS and Task 2-Microbleeds and not practically useful results yet for Task 3-Lacunes. It also highlighted the performance inconsistency across cases that may deter use at an individual level, while still proving useful at a population level.

Denoising diffusion-based MRI to CT image translation enables automated spinal segmentation.

BACKGROUND: Automated segmentation of spinal magnetic resonance imaging (MRI) plays a vital role both scientifically and clinically. However, accurately delineating posterior spine structures is challenging. METHODS: This retrospective study, approved by the ethical committee, involved translating T1-weighted and T2-weighted images into computed tomography (CT) images in a total of 263 pairs of CT/MR series. Landmark-based registration was performed to align image pairs. We compared two-dimensional (2D) paired - Pix2Pix, denoising diffusion implicit models (DDIM) image mode, DDIM noise mode - and unpaired (SynDiff, contrastive unpaired translation) image-to-image translation using "peak signal-to-noise ratio" as quality measure. A publicly available segmentation network segmented the synthesized CT datasets, and Dice similarity coefficients (DSC) were evaluated on in-house test sets and the "MRSpineSeg Challenge" volumes. The 2D findings were extended to three-dimensional (3D) Pix2Pix and DDIM. RESULTS: 2D paired methods and SynDiff exhibited similar translation performance and DCS on paired data. DDIM image mode achieved the highest image quality. SynDiff, Pix2Pix, and DDIM image mode demonstrated similar DSC (0.77). For craniocaudal axis rotations, at least two landmarks per vertebra were required for registration. The 3D translation outperformed the 2D approach, resulting in improved DSC (0.80) and anatomically accurate segmentations with higher spatial resolution than that of the original MRI series. CONCLUSIONS: Two landmarks per vertebra registration enabled paired image-to-image translation from MRI to CT and outperformed all unpaired approaches. The 3D techniques provided anatomically correct segmentations, avoiding underprediction of small structures like the spinous process. RELEVANCE STATEMENT: This study addresses the unresolved issue of translating spinal MRI to CT, making CT-based tools usable for MRI data. It generates whole spine segmentation, previously unavailable in MRI, a prerequisite for biomechanical modeling and feature extraction for clinical applications. KEY POINTS: • Unpaired image translation lacks in converting spine MRI to CT effectively. • Paired translation needs registration with two landmarks per vertebra at least. • Paired image-to-image enables segmentation transfer to other domains. • 3D translation enables super resolution from MRI to CT. • 3D translation prevents underprediction of small structures.

Benefits of a mosaic approach for assessing cortical atrophy in individual multiple sclerosis patients.

OBJECTIVE: Cortical gray matter (GM) atrophy plays a central role in multiple sclerosis (MS) pathology. However, it is not commonly assessed in clinical routine partly because a number of methodological problems hamper the development of a robust biomarker to quantify GM atrophy. In previous work, we have demonstrated the clinical utility of the "mosaic approach" (MAP) to assess individual GM atrophy in the motor neuron disease spectrum and frontotemporal dementia. In this study, we investigated the clinical utility of MAP in MS, comparing this novel biomarker to existing methods for computing GM atrophy in single patients. We contrasted the strategies based on correlations with established biomarkers reflecting MS disease burden. METHODS: We analyzed T1-weighted MPRAGE magnetic resonance imaging data from 465 relapsing-remitting MS patients and 89 healthy controls. We inspected how variations of existing strategies to estimate individual GM atrophy ("standard approaches") as well as variations of MAP (i.e., different parcellation schemes) impact downstream analysis results, both on a group and an individual level. We interpreted individual cortical disease burden as single metric reflecting the fraction of significantly atrophic data points with respect to the control group. In addition, we evaluated the correlations to lesion volume (LV) and Expanded Disability Status Scale (EDSS). RESULTS: We found that the MAP method yielded highest correlations with both LV and EDSS as compared to all other strategies. Although the parcellation resolution played a minor role in terms of absolute correlations with clinical variables, higher resolutions provided more clearly defined statistical brain maps which may facilitate clinical interpretability. CONCLUSION: This study provides evidence that MAP yields high potential for a clinically relevant biomarker in MS, outperforming existing methods to compute cortical disease burden in single patients. Of note, MAP outputs brain maps illustrating individual cortical disease burden which can be directly interpreted in daily clinical routine.